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Abstract 

The t-t' Hubbard model for the Fermi level near the Van Hove singular- 
ity is considered within the renormalization group and parquet approaches. 
The interplay of ferromagnetic, antiferromagnetic, and superconducting chan- 
nels is investigated, and the phase diagram of the model is constructed. In 
comparison with previous approaches, the account of ferromagnetic fluctua- 
tions suppresses superconducting pairing and, vice versa, the influence of the 
Cooper channel decreases the Curie temperature, so that the Stoner criterion 
is inapplicable even qualitatively. 
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I. INTRODUCTION 



Magnetic mechanisms of high-temperature superconductivity (HTSC) have become the 
subject of intensive investigations during last decades (see, e.g., Refs. It was ar- 

gued that the superconducting properties of HTSC materials are intimately related to their 
magnetic properties in the normal phase. In particular, many features of HTSC compounds 
were explained from the point of view of the competition of antiferromagnetic and supercon- 
ducting order parameters A similar physical situation takes place in ruthenates like 
Sr2Ru04, where the interplay of ferromagnetism and p-wave superconductivity is of crucial 
importance [§. Both copper-oxide systems and Sr2Ru04 are layered compounds. There- 
fore a general problem can be formulated as the investigation of the competition between 
magnetic and superconducting instabilities in two-dimensional (2D) electron systems. 

On the other hand, the problem should be concretized. There are some evidences from 
both electron structure calculations and experimental data that the Fermi surface (FS) of 
HTSC compounds at optimal doping or optimal pressure is close to the Van Hove (VH) 
singularities of electron spectrum (see, e.g., Refs. ||10|-p!^); the situation in the ruthenates 



is similar [Tj|. Due to the presence of VH singularities, the density of states at the Fermi 
level becomes logarithmically divergent, which makes the Fermi liquid unstable with respect 
to magnetic ordering or superconductivity. One may expect a priori that near VH band 
fillings the physics is determined by VH points and is not sensitive to the form of the whole 
FS. This should be correct provided that FS is not nested, since in the nesting situation 
il5| , [T6|| the contribution of fiat FS parts is also important jlTipiSl . As it was first discussed 
by Dzyaloshinskii |T^, the situation in 2D fermion system near VH fillings is similar in 
some respects to that in one-dimensional (ID) systems [^,0 or in 2D systems under the 



nesting condition p2[. It turned out that in the instability region the "normal" state can 
demonstrate a non-Fermi liquid behavior . 



The simplest model which gives a possibility to investigate the effects of VH singular- 
ities on magnetic ordering and superconductivity is the t-t' Hubbard model which takes 
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into account nearest-neighbor and next-nearest-neighbor hopping. It is often discussed also 
in connection with HTSC compounds where the value t' /t = —0.15 was determined for 
La2Cu04 and the value t' /t = —0.30 for the Bi2212 system ||2^ (we neglect third-neighbors 
hopping t" which does not lead to any qualitative changes). 

The interplay of antiferromagnetism and superconductivity near VH fillings within the 
Hubbard model with t' = was first investigated in Refs. |1T9| , |25[| . It was shown that the 
leading instability in this case is antiferromagnetic one, and the transition temperature is 



close to its mean-field value. The authors of Refs. |19| , |25[| did not include the contribution 
of the particle-hole scattering at small momenta, as well as particle-particle scattering at 
momenta near Q = (vr, vr) by the reason that these contributions are logarithmic rather than 
double-logarithmic. As we will argue in the present paper, even for weak-to-intermediate 
coupling regime these contributions should be also taken into account, which leads to an 
essential change of the results. 

Recently, the authors of Ref. have performed the renormalization group (RG) anal- 
ysis of the states close to FS (not only near VH singularities) within an approach which is 



similar to that of Ref. At t' /t = —0.30 they found competing antiferromagnetism and 



superconductivity, depending on the band filling. However, the approach of Ref. |]2^ also 
does not enable one to treat particle-hole scattering at small momenta on the equal footing 



with other contributions. As it was shown by the RG analysis in Ref. [^, the account of 
particle- hole scattering leads to occurrence of ferromagnetic phase at large enough \t'\/t; the 
criterion t' /t < —0.27 for the stability of ferromagnetism was obtained. At —0.27 < t' /t < 
it was also found that either antiferromagnetism or superconductivity takes place. How- 



ever, unlike Refs. [0,^,^, the contributions of the Cooper channel were not taken into 



account in Ref. ||28|. The backward influence of the Cooper channel on magnetic ordering 
was investigated within the T-matrix approach It was found for the non-degenerate 
Hubbard model that the Cooper channel strongly suppresses strongly the tendency to fer- 
romagnetism, so that it is possible only for t' /t < —0.35. Numerical calculations predict 
much larger values {t'/t)c for the stability of ferromagnetism at VH filling: {t'/t)c > —0.47. 



Summarizing all these approaches, one can see that we need to consider on the equal 
footing all four types of scattering to obtain the correct phase diagram, i.e. the particle- 
particle and particle-hole channels at both small momenta q and q ~ Q. First step in this 



direction was made in Refs. |31,32] within the so-called two-patch approach. The authors 



of Refs. I^T],^ wrote down approximate equations, which were very similar to the parquet 



equations in one dimension ^3|, and obtained reasonable physical results. However, they 
also neglected particle-hole scattering at small momenta and particle-particle scattering at 
q ~ Q at final stage. 

Besides that, these ID-like equations do not reproduce correctly all the peculiarities 
of the 2D dispersion law even close to VH singularities. From this point of view, most 



straightforward is the parquet approach of Refs. |ll9|-[2^ (for a review see also Ref. [0). 
It was applied to the VH singularity problem in Refs. but only the case t' /t = 

was considered (strictly speaking, this case requires an account of the whole FS because of 
nesting |r^] ). 

In the present paper we consider different phases of the t-t' Hubbard model {t' /t < 0) 
and construct the phase diagram near VH fillings within the approach |2^. Note that this 
approach is somewhat different from that used in later papers and, as we will argue 

in this work, is more correct. 

The outline of the paper is the following. In Sect. 2 we discuss noninteracting susceptibil- 
ities and consider the random-phase approximation (RPA). In Sect. 3 we consider two-patch 
equations with all four channels of scattering included and discuss the results of numerical 
solution of these equations. In Sect. 4 we consider a full parquet approach to VH problem 
and compare the results with those of two-patch approach. In conclusion we summarize 
main results of the paper and discuss possible directions of further investigations. 

II. THE MODEL AND RPA RESULTS 

We consider t-t' Hubbard model on the square lattice: 
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with 



k i 



(1) 



= — 2t(cos k.j. + cos ky) + At' cos cos ky + 4t' — (2) 

where /i is the chemical potential (we have picked out At' for farther convenience) and we 
have already absorbed the sign of t' into Eq. (|^) (hereafter we assume the transfer integrals 
t,t' to be positive, < t' /t < 1/2). 

The spectrum (^ contains VH singularities connected with the points A = (vr, 0), 5 = 
(0,7r). These singularities lie at the Fermi surface for the filling with /i = and arbitrary 
values of t' . For t' = the FS is nested, but the nesting is removed for t' /t > 0. 

Being expanded near the VH singularity points, the spectrum (^ takes the form 

= -2t(sin^ ipkl - cos^ Lpk^) = -2tk+k_ - fx (3a) 
= 2t(cos^ ipkl ~ sin^ ifky) = 2tk+k^ — /i (3b) 

where k^ = tt — k^, ky = n — ky, 

k± = sin (fk^ ± cos ipky 

k± = cos ipk^ ± simpky (4) 

(p is the half of the angle between asymptotes at VH singularity, 2ip = cos~^(2t'/t). 

At U = we have the following results for the susceptibilities at small q and q ~ 
Q = (tt, tt) (see Fig. la): 

A /(^k) " /(^k+q) ^0 , ^ X fr \ 

AB / (^k ) ~ /(^k+q) 1 ■ / c c c c \ (r:x^\ 

^q+Q = 2^ Ta—b = :T37 inin(2;Qe+, zc^i., e+e-)- (5b) 

k ^k ^k+q 

Here f{e) is the Fermi distribution function, ^± = min[ln(A/g-|-), ln(A//i)] 



za = ln[(l + VT^)/R], (6) 



and R = 2t'/t. The expressions for B ^ A are obtain by replacing ,^_|_ ^j^, 
where ^± = min[ln(A/g±), ln(A//i)]. The momentum dependence of Xq calculated with the 
spectra (^) is shown in Fig.2. Since both the susceptibilities are divergent, we have at least 
two competing order parameters. In fact, two other polarization bubbles of Fig. lb, which 
are responsible for zero-momentum and 7r-pairing, are also divergent at small q: 



where 



Hi 



k+q 



Co 

2ttH 



(7a) 



n- 



AB 
q+Q 



1 - fjs^) - /feq 
^ ^k+q 



2ttH 



min(^+, 



(7b) 



Co 



1/Vl cq = tan-\R/Vl- R'^)/R 



For Ilq we again have the replacements A ^ B and ^_|_ ^_|_, ^_ — ^_ in ([ra|). 

In the RPA the expressions for particle-hole and particle-particle susceptibilities read 

Xq 



l-f/Xq 

nq 

1 + UUa 



(8) 
(9) 



Thus n decreases when the Coulomb interaction is taken into account, while x increases 
and can diverge at some U. In particular, we have the conventional Stoner criterium of 
ferromagnetism Uxo = 1, or 



^^°ln^^l 



(10) 



2TrH p 

where p = max(T/t, p/t, A/t) (A ~ S* is the spin splitting). The solution to this equation 



reads 



p = Aexp \-2^r^^/t^^^^J2fy/U 



Therefore the ferromagnetism is present at any U ; moreover, at f/ ~ 2tt [t —(2t') ] ' one can 
expect that it becomes saturated. Similarly, considering the antiferromagnetic instability 
we obtain 
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u 

2^t 



A 



A. 



min(ln — , zq In — ) = 1 



(11) 



which gives 



p = A 



exp(-yM7f7), U/{2txH) > 1/4 
exp{~27rH/zQU), U/{27vH) < 

so that antiferromagnetism is favorable at small t'/t. 



However, as it was discussed first by Dzyaloshinskii and coworkers [|T9]-|2^ , RPA is incor- 
rect even in the weak-coupling limit, except for the case when only one bubble is divergent. 
Since in the VH case all four bubbles of Fig.l are divergent, we have to use the parquet 
approach ||l9|-p2| instead of RPA. While in ID case the parquet equations reduce to conven- 



tional differential RG equations (see, e.g., Refs. [0,^,^), for higher space dimensionalities 



we have coupled integral equations. First we consider the approach of Refs. ||31|j34 which 
uses mapping of the full parquet equations on an "effective" ID problem, i.e. so-called 
two-patch equations. 



III. TWO-PATCH EQUATIONS 

The authors of Refs. [^,^ proposed the approach which neglects the difference between 
^+ and (and consequently between ^+ and ,^_) and introduced a single scaling variable 
^ = min(^+, ^_). Note that this approach is not strict, in particular because of the 

presence of double- logarithmic terms in (|5bD and (^). At the same time, as we will see below 
(see also Ref. [^), this reproduces correctly main features of the exact parquet equations. 

The two-patch equations read 

7^ = 2c?i(07i(72 - 7i) + 2^27174 - 2 ^37x72 

72 = MOill + 73) + 2^2(71 - 72)74 - dsill + 72) 

7^ = -2rfo(07374 + 2rfi (073(272 - 7i) 

7; = -^^o(0(73 + 74) + ^^2(7? + 27172 - 272' + 74) (12) 
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where 7- = d'ji/d^, 



do{0 = 2co^; diiO = 2mm{^, zq) 

d2 = 2zo; dz = 2cq (13) 

and four vertices 71-4 are defined in Fig.3. In these notations, 7i(0) = go = U/{4:7iH) 
corresponds to the Hubbard model. While only the case ^2,(^3 do, di was considered 
in Refs. [^,^, we perform a more general consideration where all the bubbles are taken 
into account. We have also taken into account the coefficient cq to treat correctly the t' 
dependence of the amplitude of particle-particle scattering. Note that the equations (|T2p 
are very similar to those in the ID case with the difference that in the latter case 

one has 



do = d2 = 0; di = ds 



(14) 



The complete discussion of the physics of the equations (12) in the ID case is given in 



Ref. ||2T[. In two dimensions, the coefficients do and di become ^-dependent because of 
the presence of double-logarithmic terms. As we have already mentioned, this gives only 
approximate treatment of such terms. The equations ([12D give a possibility to investigate 
the interplay of all the four scattering channels. 

For ferromagnetic, antiferromagnetic, and d-wave superconducting susceptibilities we 
have 



XF,AF,<i-Sc(0 = / '^C'^2,l,o(C)^^,AF,d-Sc(C) 



(15) 



where T satisfies the equation 

d In Tp A F,d-sc 



c^2(7i + 74) 

< ^1(0(72 + 73) \ (16) 
c?o(0(73 -74) 

V. 

Thus, when 71 and 74 are simultaneously relevant, we have ferromagnetic ordering, while 
72 and 73 lead to antiferromagnetic ordering. For the superconductivity, we have a more 



complicated combination of relevant and irrelevant vertices. 



The results of the solution of equations (p!2D-(p^) for various values of go and t' / 1 are shown 
in Fig. 4. Depending on the values of go and t' /t, ferro- or antiferromagnetic susceptibility, or 
(i-wave superconducting response diverges first. The parameter dependences of the critical 
energy scale fic = texp{—C,c) are shown in Fig.5. This scale can be approximately identified 
with the critical chemical potential or transition temperature. For comparison, the RPA 
results for the stability of ferro- and antiferromagnetism are shown too. One can see that 
the values of transition temperatures obtained from the two-patch equations are much lower 
than the corresponding RPA results. 

To understand qualitatively the nature of the critical temperature lowering, we may 
neglect the interpatch scattering. In the ferromagnetic case we have only one nonzero vertex 
74, and the equation for it has the form 

7; = -2{zo - CoOll (17) 

so that 

= 1 + g,{coe - 2zoO ^^^^ 

The modification of the Stoner criterion takes the form 

goi2zo\n- - coln^ -) = 1 (19) 
Then we have from (|19D at Zq = Cq (which is the case of t-t' model) 



A ^/zogo^ 1 
m — = 1 > 



IJ, zogo 2zogo 

(for R close to unity we have zogo > !)• Thus, the decrease of the Curie temperature in 
comparison with the mean-field approach (cq = 0) is directly connected with the account 



of the Cooper bubble, which is in agreement with the T-matrix approximation Note, 
however, that the structure of Eq. (0) is different from that obtained in the T-matrix 
approach. 

The resulting phase diagram in U-t'/t plane with all the scattering channels being in- 
cluded is shown in Fig. 6. One can see that the cZ-wave superconducting response is strongly 
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suppressed by particle-hole scattering processes. This is the leading divergent response only 
at small values of coupling constant qq < 0.04, which corresponds to U < 1.6t. The critical 
temperature in this region is also exponentially small (Tc ^ fic ^ 6xp(— 1/(70) )• 



IV. PARQUET EQUATIONS 

Now we pass to the consideration of the full parquet equations and compare the results 
of their solution with approximate RG equations of Sect. 2. We use the generalization of 
the approach of Ref. to the case of two dimensions. 

In the parquet approach (see Appendix) we have for each vertex i = 1...4 three types of 
bricks, which are shown in Fig. 7: the Cooper brick Ci{^±,ri±), and two zero-sound bricks, 
Zi{^±, rj±) and Zi{^±, r]±). Up to a logarithmic accuracy, they depend on ^± = \n{A/k±) and 
r]± = ln(A/g±) only, k± = ki± + k2± and q± = maxjfcs-i- — ki±, k^± — k2±} being the Cooper 
and zero-sound momenta transfer. The vertices 7i(^±,''7±) in different regions of ^± and ri± 
are given by 

7i(^±, Tl±) = li{^±^ V±) =90 + Ci{^±, r]±) + Z,{r]±, r]±) + Zi{r]±, r]±) > r]±) 

7i(e±, V±) = 7-(e±, V±) ^90 + a(e±, e±) + r^i) + Z,(e±, ^i) < V±, V± < V±) 

7i(e±, v±) = 7-(e±, v±) ^90 + m±, e±) + ^.(e±, e±) + m±, v±) (e± < v±, > ^1'^) (20) 

where r/^'^'' = ln(A/|fc3± — ki^2±\)- Following to Ref. [^, we have taken into account that at 
^± > V± the Cooper brick depends on both ^± and ?7±, while the zero-sound bricks depend 
only on ri±. Vice versa, at C,± < f]± the Cooper brick and one of the zero-sound bricks depend 
only on ,^-|-, and another zero- sound brick depends on both ^± and ?7±. 
When all momenta are of the same order, i.e. = ti±, the vertices 

7.(e±,e±) = 7.(e±) (21) 

are analogous to those introduced in Sect. 3 with the only difference that now they depend on 
two scaling variables However, unlike the ID case, the parquet equations do not reduce 
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to the equations for 7j(^±), but contain the full dependence 7j(^±,?7±). The corresponding 
equations are presented in Appendix. As discussed in Appendix, the approach we use gives a 
possibility to treat the 2D situation in a more correct way in comparison with the approach 
of Refs. 

The parquet equations were solved numerically. To this end, we placed the variables ^, rj 
on a grid with 16 points in each dimension, so that the total number of vertices to be taken 
into account is 3 ■ 4 ■ 16^ ~ 8 ■ 10^. It is important that the grid was chosen for the logarithmic 
variables ^, i], but not for the momenta themselves. This gives a possibility to use simple 
integration methods (e.g., the trapezium method) to obtain the results which are correct to 
logarithmic accuracy. The resulting system of 8 ■ 10^ algebraic equations was solved by the 
Zeidel method. 

The structure of the solutions of the parquet equations is quite similar to that in the 
two-patch approach, except for that now we have momenta-dependent vertices. Again, 
the relevance of 71 and 74 with ^± = ''7± = ^ leads to ferromagnetic ordering, while the 
relevance of 72 and 73 at ^± = r^-i- = ^ to antiferromagnetic one. The results of solution 
of the parquet equations are shown in Figs. 8, 9. One can see that the results coincide 
quahtatively with those of the two-patch parquet approach of Sect. 3. At not too large t'/t, 
the antiferromagnetic instability occurs first, while for t'/t close to 1/2 the leading instability 
is ferromagnetic one. The superconductivity occurs also only for very small Qq. 

The transition temperatures obtained within the parquet approach are larger than those 
obtained from two-patch equations, but are still lower than the RPA results. In particular, 
in the limit of small t'/t the parquet calculations [|19|,^, which do not take into account 
single-logarithmic contributions of the loops Fig.la,d, yield for go = 0.1 the critical value for 
stability of antiferromagnetism = 5.2, which is close to RPA result, ^"^ = 5.0. At the same 
time, our parquet calculations give larger value, = 6.37 (the result of two-patch equations 
is = 18.2). The region of stability of d-wave superconducting phase is even smaller than 
that obtained from two-patch equations. 

The critical concentrations ric for the stability of ferro- or antiferromagnetic and super- 
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conducting phases close to VH filling can be estimated from the critical chemical potential 
with the use of the condition 



n 



E/(^i') (22) 



Using the form of spectrum @ and taking the limit of filling close to VH one, we obtain for 
< t 

e ^ic , At ^cexp(-^c) .r>o\ 

oric = ric — nvH = , m ■; — - ~ , (zS) 

where rivH is the VH filling. In particular, for qq = 0.1 {U = 3.95t) we have from Fig. 8 

5n^ = 0.01 (AF phase, t' /t 0) 

dric = 0.03 (F phase, t'/t = 0.45) (24) 

Thus, except for the limit R ^ 1 {f t/2), the critical concentrations are very small, which 
is in qualitative agreement with the results of Ref . [E^] . Because of the exponential smallness 



of the critical chemical potential, the critical concentrations for the superconducting phase 
are even smaller than those for the magnetically ordered phases. 



V. CONCLUSION 

Now we summarize the main results of the paper. Using the two-patch equations (Sect. 
3) and parquet equations (Sect. 4) we constructed the phase diagrams of t-t' Hubbard model 
(Figs. 5,6,8,9) at the fillings which are close to Van Hove one. It was argued that the simulta- 
neous account of all the scattering channels is important in considering the VH problem, the 
smallness of contributions of some channels (logarithmical vs. double-logarithmical diver- 
gence) being compensated by the growth of relevant couplings. Both the approaches used, 
two-patch and parquet ones, give similar phase diagrams. In agreement with the previous 
approaches p8|^0|, antiferromagnetism is favorable for small t'/t, while ferromagnetism for 
larger values of t'/t. The stability of antiferro- and especially ferromagnetism is greatly re- 
duced in comparison with the corresponding mean-field criteria. Thus the Stoner criterion is 
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completely inapplicable for the systems with VH singularities; depending on the value t' /t, 
it overestimates the critical temperature by 2-10 times. This conclusion is in qualitative 
agreement with the results of Ref. |^. Besides that, the mean-field approach is unable to 
determine the critical value {f /t)c which separates the ferro- and antiferromagnetic phases. 

Unlike Ref. (t'/t)c is [/-dependent and decreases with increasing U. Although the RG 
(and also parquet) approach is unable to describe the ordered states, from scaling arguments 
we have S oc (fic/t)^, where P is the magnetization critical exponent. With increasing U, 
the ferro- and antiferromagnetic states are characterized by large magnetic moments, and 
ferromagnetism possibly becomes saturated. However, these values of U are not described 
by perturbative approaches and should be treated in the strong-coupling hmit. At the same 
time, determining parameters of the ferro-antiferromagnetic quantum phase transition would 
be of interest, especially the critical exponents. One can expect that they are independent 
of the coupling. 

Another result of the paper is that the tendency to d-wave superconducting pairing is 
considerably reduced in comparison with the treatments of Refs. [^,^ : it can occur only at 
very small values of U. Of course, this concerns only the pairing due to the VH singularities 
themselves; the pairing can be further enhanced by other factors. This can be also the 
subject for future investigations. Details of the electron spectrum, especially the form of the 
fermion Green's functions close to the phase transition into the ferro- or antiferromagnetic 
state are beyond the scope of the present paper. Although the marginal and non- Fermi- 
liquid behavior p3[ was found (see also the discussion in Ref. ||3^), this problem needs 



further investigations, since a simultaneous account of all the scattering channels can be 
important in this case too. For example, only one of four scattering channels was included 
in Ref. |2| . 



We believe that the results of the present paper can be also important for the theory 
of itinerant-electron ferromagnetism. A standard consideration (including contemporary 
spin-fluctuation theories |^) starts from the RPA approach. It was noted in Ref. 



that for almost all known itinerant electron ferromagnets the Fermi level lies near a 2D-like 

13 



VH singularity. This is a result of merging two weaker 3D square-root singularities along 



symmetrical directions in the Brillouin zone . We have shown that under such conditions 
the RPA approach and the Stoner criterion are not applicable even qualitatively because of 
the strong interference with the Cooper channel. Of course, the effect of the logarithmic VH 
singularity in the 3D case is not exactly the same as in the pure 2D case considered here, so 
that the 3D problem needs further investigations. However, the naive Stoner criterion is in 
any case doubtful and needs a careful justification. 

In this respect, it would be interesting to generalize the results of the present paper (at 
least those from the two-patch equations) on the degenerate-band Hubbard model. As was 
argued in Ref. [^, in this case the suppression of ferromagnetic ordering is much weaker 



than for the nondegenerate model considered. One can also expect that the particle-hole 
scattering with small momenta will not renormalize superconducting channel as strongly as 
for the non-degenerate model. However, these statements need further justification since 
the diagram series in the degenerate and non-degenerate cases look like rather similar. 

The research described was supported in part by Grant No. 00-15-96544 from the Russian 
Basic Research Foundation (Support of Scientific Schools). 

APPENDIX A: THE PARQUET APPROACH IN ONE AND TWO DIMENSIONS 

First we consider the simple model of spinless fermions 

H = ^kcick + ^(p)cick_pcick+p (Al) 

k k,p 



with gikp) = go- In one dimension we have the representation for the renormalized 
vertex 7 (Fig. 5) 

7(e, V) = 9^ + ri) + Zi(^, r/) + Z^ii, r,) (A2) 
where the bricks are given by 

C(e,r/) = -c /" rfC7'(C:^)7'(C,^), e > V 





14 



Zi,2{^,v) = Zi,29o j c?C7i,2(C:e)7l,2(C,^), e < r] (A3) 



and ^ = ln(A/|fci + ^2!); f] = ln(A/ maxjfca — ki, — ^2})- Here 

and we assume that the Cooper, ZS and ZS' loops are logarithmically divergent with the 
coefficients c and 2:1,2 respectively (we generalize here the approach of Ref. ||2^ to the case 
where both the channels, ZS and ZS' contain divergences). The vertices in ([A3| ) are given 
by 

7'=(0=^?o + ^i(e,0 + ^2(e,0 

7M(0 = ^?o + c^(e,o + ^2,i(e,o (A5) 
7'(e,r/) = 7^(0 + c^(e,r/) 
75 2(e,r/) = 71,2(0 + ^1,2(^,^7) 



The equations (|A2|) - ([A5|) form the closed system of parquet equations for the ID spinless 
case. The validity of these equations can be demonstrated for the trivial case 2:1 = 2:2 = 
where the direct ladder (RPA) summation is possible. In this case Zi^2{^, v) = 0; 7*^(0 = 9o, 
and we obtain from (|A^), ( |A5| ) the standard ladder equation 

7(^, V) = 9o- c^o^7(^, V) (A6) 

This has the solution 

7(e,r/) = -^ (A7) 



Now we return to the general case zi, ^2? c 7^ 0. As it is shown in Ref. [^, the equations 
-(|A5D at ^ = ?7 can be reduced to 

7(e, = 7(0 =90 + {z, + Z2- Co) j l\C)dC (A8) 
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which is equivalent to the differential RG equation 

% + Z2 - c)7' (A9) 



Since in one dimension zi = c, Z2 = 0, 7 is marginal and we have a Luttinger-liquid behavior 
p7| , |33|jr7| . Alternatively, the equations ( [A8| ) or (|A9| ) can be obtained directly with the use 
of the Sudakov's trick |^0| or standard RG approach |^ without considering the general 



dependence 7(^,^7). Thus, in the ID case the parquet and RG approaches are equivalent. 

In the 2D case we have two pairs of variables, ^± and ri±. However, unlike the ID case, 
now two possibilities occur: the momentum integration in bubbles can be logarithmical 
or double-logarithmical. For example we consider the case where the integration in the 
Cooper bubble is double-logarithmic while in the zero-sound channel this yields only simple 
logarithms (which is similar to the situation for VH singularities). Then it can be checked 
by a direct comparison with perturbation theory that the equations 



C{i±,r]±) = -ego J j rfC+ciC-7'(C±,^±)7 (C±;^±), C± > V± 


Zi,2i^±,v±) = zi,2go j rfC+7i,2(C+,'^+,^-)7l,2(C+,^-;^±) 


V- 

+zi,2go j rfC-7i,2(^+> C-, ^-)7l,2(^+> C-; v±), ^±<v± (Aio) 



with 



i\i±.v±)=i\v±) + c{i^,v±) 

l\{i±,V±) = lt{i±) + ^.(e±,r/±) (All) 

give the parquet solution of the problem. Note that beyond one dimension the integral 
parquet equations do not reduce to differential ones. 
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The above approach is different from the approach of Refs. where a standard 

RG scheme was apphed in one dimension, while momentum dependence in another dimen- 
sion was taken into account exactly rather than to logarithmic accuracy. However, the 
applicability of last approach is doubtful. Indeed, in the ID case the equation ( |A8| ) can be 
considered as a logarithmic approximation to the Bethe-Salpeter equations 

A 

C{ki, /c2, h) = go- cqQq / dk'-i{ki + k,k2- k, ks) 





A 



Zi{ki, k2, /cs) = go + ziQo J dk'j{ki, k-^ + /c, k^) 



A 

^2(^1, ^2, fcs) =90 + Z2go j dk-f{ki, k2 + fc, /cs + k) (A12) 







If we have both slow and fast momenta we need to combine (^^) and ( |A12| ) which is impos- 



sible since ( [A 81 ) is quadratic in 7 while ( |A12| ) is linear. The equations of Refs. [Q|T9|,|T7 

? A 

C{ki, k2, ks, = fi'o - Co j dQ j dk'y{ki, k2, k^ + k; C)7(^i + fc, ^2 - k, ks; Q 



? A 

Zi{ki,k2,k3,^) = go + zi j dC j dk-f{ki, k3 + k, /cg; C)7(^i + ^, k2,k3 + k;C) 



5 A 

Z2{ki, k2, fcs, = 9o + Z2 j dC j dk'j{ki, k2 + k,k3 + k; Qlih + k, k2, ks; () (A13) 



are not fully correct. If we suppose that 7 does not depend on ^, we do not reproduce the 
ID Bethe-Salpeter equations ( |A12| ). At the same time, the approach of Ref. I^O] is free from 
these problems. 



The generalization of equations ( |A10| ) to the full VH problem is trivial. The parquet 
equations have the form 



Ci(e±; v±) = -CQ / [7?(C±r7±)72'(e±; C±) + 72(C±^±)7i (e±; C±) 
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+73lC±e±)7^(e±; C±) + 7f (C±e±)7^(e±; C±) - 27|(C±e±)7^(e±; C±))] 

^i(e±; V±) = Zo [ 7l(C^)7l (C±, V±) + Z0 [ 7i'(C^)7l(C±, V±) 

C2(e±; r^i) = -CQ / K(C^)7f (e±, C±) + 72(C±^)72 (e±, C±)] 

^2(e±; V±) = zo r [7l(C^)72(C±, V±) + 7l(C±^)7l(C±, r/±) - 27l(C±^)7^(C±, r?±)] 

+^0 b!ic^hiic±,v±) +iiic^Hic±,v±) - 272ic±e^)7i(c±,^±)] 

Jzo 

m±; V±) [ [7l(C±e^)7^(C±, + 7f (C^)7^(C±, v±)] 

J ZQ 

Cs{^±;v±)--co r74(C±;^)73(e±;C±)-co ^ Y3(C^H{C±;C±) 

J CO J CO 

Zs{^±; V±) [ [73(ae^)(72(C±, v±) - lt{C±, v±)) + hliC^) - Yi(C^)hi{C±, v±) 

J ZQ 

+ ll{C^){li{C±,V±) - liiC±,V±)) + illiC^) - 73(C^))7l(C±,r?±)] 

m±; v±) -ZQ [ [i!(c^mc±, v±) + 7l(C^)7i(C±, v±)] 

J CO J CO 

Z,{i±\ V±) = zo r bl(C^H{C±: V±) + 7l(ae^)7i(C±, V±) - 27l(ae^)7i(C±, V±)] 

J zo 

+Zo f [7l(ae^)7l(C±, V±) + 7f (C^)72(C±, V±) - 272 (C^)72(C±, V±)\ 

J Zo 

Z,{^±;V±) = Zo r7l(C^H{C±:V±) + Zo f ll(^)l{{C±.ri±) (A14) 

J zo J zo 



The vertices are now given by 



li{v±) = 9o + Zi{r]±, r/±) + Zi(r/±, r]^) 

7f (e±) = 90 + a{^±, e±) + ^^(e±, e±) 

7^(e±, V±) = 7f (e±) + m±, V±) (A15) 
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and the regions of integration are defined by 

^ /(C+,C-)= I J rfC+c?C-/(C+, C-) 



zo 



where 



^ /(C+,C-)= J I rfC+c^C-/(C+,C- 



15+1 -15-1 

[5+1 / max{0,d''} max{-15_|signC+,ci^>}\ 



-15+1 VminiOX^^H min-f-lf_lsien 



\mm{0,d'^} min{-|e-|signC+,d''} / 



dC+dC- 

COS 2(/7 



max{0,|5+|sign5_} 



+ / dC 

min{0,|C+lsign5_} 
max{0,|5-|sig<+} 

+ / * 

mm{0,|5-|sigii5+} 
W+l l»7-l 



min(cos29?/c+,p_) +]?_ 



cos 2(pk+ + sin 2(/7p_ 
min(cos2(/7A;_,p+) +p+ 



cos 2(/7A;_ + sin 2(pp+ 

I 

£f {(+,(-)= I dC+f{C+,V-)+ I dC+f{r,+,C-) 



kl + kl + 2 cos 2(fik+k. 

/(C+,e-) 

/(e+,c-) 



/(C+,C- 



-\v+\ 



\v-\ 
\v-\ 



/(C+,C-)= dC+fiC+,v-)+ dc+fiv+x-) 



-\v+\ 



^ max{0,|»7_|signC+} max{ci^\ci^' }^ 



-\v+\ 



min{0,|7j_|signC+} 



dC+dC- 



kl + ki + 2 cos 2(fik+k- 



f {(+,(-) (A16) 



fc± = Asign(C±) exp(-|C±|) 
L{k±) = sign(A;±) In |A/A;±| 
= -L[k+/ cos 2ip] 
-L[k^/ cos2ip] 
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Figure captions 

Fig.l. Diagrams (bubbles) for noninteracting susceptibilities near q = and q = Q in 
(a) Peierls channel (b) Cooper channel. Solid and dashed lines correspond to the electron 
Green functions near A and B singularities with the spectra ( pa]) and (pED respectively. 

Fig. 2. The momentum dependence for (a) noninteracting susceptibility and (b) nonin- 
teracting Cooper response at t'/t = 0.3 

Fig. 3. The vertices ji [i = 1...4). The solid lines inside the circles show which incoming 
and outgoing particles have the same spin projection. 

Fig. 4. The ferromagnetic (solid line), antiferromagnetic (long-dashed line), and d-wave 
superconducting (short-dashed line) susceptibilities for the two-patch model with (a) t'/t = 
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0.15; qq = 0.10 (b) t'/t = 0.45; go = 0.10 (c) t' /t = 0.30; go = 0.01 

Fig. 5 The phase diagram for the two-patch model in ^-t' /t coordinates for g^, = 0.1 
{U = 3.95t). Dotted hne is the mean-field boundary for antiferromagnetic phase, dot-dashed 
line for ferromagnetic one. 

Fig. 6. The phase diagram for the two-patch model in go-t'/t coordinates at Van Hove 
filling (/x = 0). 

Fig. 7. The representation of the vertex in the parquet approach as a sum of bricks for 
the Cooper (C) and zero-sound {Z, Z) channels. 

Fig. 8. The phase diagram from parquet equations in /i-t'/t coordinates for go = 0.1 
{U — 3.95t). The lines are the same as in Fig.5. 

Fig. 9. The phase diagram from parquet equations for Van Hove fiUing (// = 0) in go-t'/t 
coordinates. 
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This figure "fig2a.gif" is available in "gif" format from: 
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This figure "fig2b.gif" is available in "gif" format from: 
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